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I.  INTRODUCTION 


In  1946  V.  L.  Ginsburg  and  I.  M.  Frank  predicted  the  existence  of  transition 
radiation  [Ref.  1:  pp.353-362].  This  effect  was  discovered  during  a  period  in  which 
Russian  theoreticians  made  an  extensive  effort  in  the  area  of  radiation  producing 
phenomena  after  the  discover}'  of  the  Cerenkov  effect.  Actually,  transition  radiation 
can  be  considered,  from  a  phenomenological  point  of  view,  as  a  Cerenkov  effect  of  the 
seond  order  since  it  occurs  when  a  charged  particle  crosses  the  boundary  between  two 
media  with  differing  dielectric  constants.  More  generally,  this  effect  will  take  place  in 
the  presence  of  the  inhomogeneities  in  the  medium. 

Through  the  years,  experimental  investigation  devoted  to  this  effect  were 
undertaken  primarily  with  nonrelativistic  electrons,  in  order  to  firmly  establish  the 
validity  of  the  theory,  in  spite  of  inconvenience  arising  at  these  low  energies  from  the 
simultaneous  existence  of  scattering  of  the  electrons  and  of  bremsstrahlung  radiation. 
Interest  in  this  type  of  radiation  was  renewed  after  G.  M.  Garibyan  predicted  that  the 
transition  radiation  yield  was  propotional  to  the  energy  of  the  panicle.  The 
application  of  this  propeny  to  the  detection  and  identification  of  individual  panicles 
was  the  main  goal  of  high  energy  physics  instrumentation,  cosmic  rays,  and 
astrophysics. 

In  recent  years  the  in-phase  addition  of  radiation  from  multiple  boundaries  has 
been  predicted  and  observed  [Ref.  2:  p.3594).  For  a  single  interface,  the  radiation  yield 
produced  by  an  individual  electron  is  very  weak;  of  the  order  of  the  fine-structure 
constant  /.  e.,  roughly  one  photon  for  a  hundred  electrons.  For  M  boundaries  the 
spectral  intensity  varies  as  M  ,  so  that  the  emission  can  be  much  greater  than  for  a 
single  interface.  In  this  thesis  the  effects  of  light  absorption,  and  random  variation  in 
foil  thickness  are  included  in  the  design  of  transition  radiation  sources  for  the 
generation  of  x-rays.  The  foil  dimensions  must  be  held  to  close  tolerance  in  order  to 
satisfy  the  coherence  conditions  from  each  boundary. 


II.  THEORY 


A.  IN-PHASE  ADDITION  OF  THE  TRANSITION  RADIATION 

Transition  radiation  occurs  when  a  moving  charged  particle  encounters  a  sudden 
change  in  dielectric  constant  at  the  interface  between  different  media  (  e.g.,  between  a 
vacuum  and  a  solid).  Ordinarily,  a  particle  which  moves  with  constant  velocity  does 
not  radiate  unless  it  is  in  a  refractive  medium  and  the  particle  exceeds  the  velocity  of 
the  electromagnetic  wave  in  medium  --  as  is  the  case  for  Cerenkov  radiation.  However, 
if  the  interaction  length  is  limited  ,  or  equivalently,  if  the  dielectric  constant  changes 
suddenly,  then  velocity  matching  is  not  important.  The  minimum  distance  over  which 
an  electromagnetic  wave  and  a  charged  particle  can  exchange  energy  is  called  the 
formation  length  and  is  given  by  [Ref.  3:  pp.  40-44] 

2cP 

Z.  -  - ,  (2.1) 

co(l-pv'Vsm-ff) 

where  y  =  (l-P2)*'2’  £;  (i=l,  2)  =  1  -  (o.  to)2  are  the  dielectric  constants  of  the  two 
media,  oy  are  their  respective  plasma  frequencies,  P  »  v  c,  v  is  the  speed  of  the 
electron,  c  is  the  speed  of  light,  hto  is  the  photon  energy,  and  0  is  the  angle  cf 
emission.  For  relativistic  electrons  P  — ^  1  ,  sin0  0,  and  1  y2  ~  2(1  -P),  so  that 

4kp 

zi  ~  ’  <—> 

(l,y)'  +  02  +  (toi  to) 


where  "k  =  c.to. 

In  traversing  an  interface,  the  number  of  photons  per  unit  time  emitted  by  an 
electron  is  proportional  to  the  dot  product  of  the  particle  velocity  and  electric  field 
strength.  For  a  single  interface  V.  L.  Ginzburg  and  I.  M.  Frank  calculated  the 
differential  production  efficiency  for  transition  radiation  as  [Ref.  1:  p.353] 


d2N  a02o>  , 

— ^  =  [ - T T-KZ.-Z,)2, 

dl2da>  16tt‘c‘  1  * 


(2.3) 
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where  d“N  dftdo)  is  the  number  of  photons  emitted  per  frequency  per  solid  angle  by  a 
single  particle  crossing  a  single  dielectric  interface,  a  =  1,1 37  is  the  fine-structure 
constant,  ft  is  the  solid  angle  in  steradians,  and  Zj  and  Z,  are  given  by  Eq.  (2.1) . 

The  spectral  intensity  f:  )m  an  electron  crossing  M  foil  pairs  perpendicularly, 
each  composed  of  two  materials  of  thickuess  and  £,  can  be  calculated  by  considering 
the  field  amplitudes  at  some  observation  point  [Fig.  (2.1)].  The  results  should  be 


simply  the  coherent  sum  of  the  radiation  amplitudes  from  2\I  single  interfaces.  The 
phase  ractor  must  be  properly  included,  taking  into  account  that  the  particle  traverses 
different  interfaces  at  different  times. 


where  N  is  the  number  of  photons  emitted  by  a  single  particle  traveling  through  a 
dielectric  stack,  where  <r  =  Jij  ^  +  M2’  x*  VZl  +  VZ2’  and  and  f*2  are  the  linear 
absorption  coeflficents  of  the  spacing  and  foil  media,  respectively  (jij  —  0  for  vacuum 
spacing).  When  the  periods  of  F(M,X)  are  not  experimentally  resolvable  and  the 
photon  absorption  of  radiation  is  not  small,  Eq.  (2.5)  becomes  [Ref.  5:  p.918] 


<F(M,X)> 


1  -  exp(-M<r) 
<J 


(2.6) 


At  high  photon  energies,  where  x-ray  absorption  is  small  or  negligible,  so  that  (<r~0) 


F(M.X) 


sin2!  MX) 
sin2X 


(2.7) 


Then  Eq.  (2.4)  can  be  written  as 

d2N  d2N  ,  „  sin2(MX) 
dfidw  di2dto  2  2  sin"X 


(2.8) 


When  the  spacing  between  the  foils  Ej  exceeds  the  formation  length  of  the  gap 
material  Zj,  F(M,X)  varies  rapidly  compared  with  the  single-interface  term  given  by 
Eq.  (2.4)  and  the  peak  spectral  intensity  is  found  to  vary  as  the  square  of  the  number 
of  interfaces  when  X  =  nt  and  r  is  an  integer,  neglecting  photon  absorption  in  the  foils, 
variation  in  foil  thickness  and  the  multiple  scattering.  From  Eq.  (2.S)  it  is  possible  to 
determine  the  foil  thicknesses  that  maximize  photon  production.  This  occurs  for 

X  =  rtl  ,  (2.9) 

and 


C2  Z2  ~  (m  -  1/ 2)jt  , 


(2.10) 


N 
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Where  r  and  m  are  positive  integers.  Substituting  these  conditions  into  Eq.  (2.S)  gives 


d2\ 

dftdo 


4M2I- 


d2\ 


d&d<0 


■I- 


(2.11) 


Eqs.  (2.9)  and  (2.10)  are  the  coherence  conditions,  /.<?.,  the  requirements  for  in-phase 
addition  of  the  radiation  from  all  interfaces.  To  minimize  photon  absorption  and 
electron  scattering,  Cj  and  t2  should  be  as  small  as  possible  which  means  choosing  r  = 
m  =  1,  for  these  values. 


(2.12) 


Figure  2.2  Condition  for  Coherent  Phase  Summation. 

With  the  assumption  that  c.  ->  ~  1,  from  Eqs.  (2.1)  and  (2.9)  it  is  found  that  the 
angle  for  which  resonance  is  satisfied  is  given  by  (Ref.  6:  pp. 264-2^5], 


COS0 


[i  P  - 


271  rc 


12.13) 
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The  in-phase  radiation  adds  coherently  and  the  resonance  transition  radiation  is  shown 
diagramatically  in  Fig.  (2.2). 


Figure  2.3  Schematic  Diagram  of  the  Transition  Radiation  Cone. 


B.  TRANSITION  RADIATION  CONE 

Fig.  (2.3)  shows  the  transition  radiation  emitted  in  a  tight  forward  direction. 
There  is  also  an  apex  angle  that  maximizes  the  radiation  d"N  dlidto  from  a  single 
interface.  For  the  conditions  0:  <  <  1.  Cj  •>  ~  1.  and  P  —  1.  this  optimizing  angle  can  be 
derived  from  Eq.  (2.13)  [Ref.  4:  p.  1290] 

02oP,  =1  3{-(6,+52)  +[(61  +62)2+  125152]'/j}.  (2-14) 
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where  5;  =  [( 1/y)2  -H (<o.  co)2]/2  (i=  1.2) .  If  both  Eq.  (2.1)  and  Eq.  (2.13)  are  satisfied 
for  0  =  9opt,  the  d2N/dftd<0  is  optimized.  In  general,  (02>  >0^,  where  C0j  :s  the 
plasma  frequency  of  the  medium  in  which  the  foils  immersed.  For  the  case  of  a  toil 
stack  in  vacuum.  (0j  ~0  and  8opt  reduces  to 

9  ~  1,'y  .  (2.15) 

opt  1 


As  the  electron  energy  increases,  the  cone  apex  angle  becomes  small.  This  gives  a 
highly  directional  and  intensive  forward  beam. 


Figure  2.4  Coherent  and  Incoherent  Transition  Radiation  from  a  Foil  Stack. 

For  incoherent  emission  the  photons  can  be  separated  from  the  emitting  charged 
particle.  For  coherent  radiation  where  there  is  inter-foil  phase  addition  of  the 
radiation,  the  lobe  pattern  will  be  broken  into  smaller  lobes  due  to  interfoil  coherent 
addition  of  soft  x-rays  (Fig.  (2.4)]. 


C.  COMPUTER  SPECTRUM  ANALYSIS 


By  taking  the  approximate  multi-foil  absortion  [Eq.  (2.6)],  the  analytic  expression 
for  the  transition  radiation  [Eq.  (2.4)]  can  be  written  as 


d2N 

dfidti) 


d2N  ,  .  l-exp(-M<J) 

4  [i?rr-]sin{l  2;Z^ — - — i 

dfldto  1  *  <T 


(2.16) 


This  expression  was  used  for  the  differential  production  efficiency  in  spectrum  analysis 
below.  The  analytic  expression  can  be  analyzed  by  dividing  it  into  the  three  parts: 

1.  The  single  interface  intensity  Jd2X0  dCldm  [Eq.  (2.3)] 

2.  The  single  foil  interference  dependence  sin2^l ,  Z,)  [Eq.  (2.4)] 

3.  The  multi-foil  interference  dependence  F(\l,X)  [Eqs.  (2.5)  and  (2.6)] 

1.  Formation  Length 

The  key  to  understanding  the  dependence  of  the  three  main  terms  on  the 
photon  energy  and  the  stack  variables  is  understanding  the  dependence  of  the 
formation  length  [Eq.  (2.1)]  on  the  photon  energy  and  stack  variables.  The  dependence 
of  the  transition  radiation  intensity  on  the  formation  length  can  be  divided  into  three 
parts. 

a.  The  single  interface  intensity  depends  on  (Zj  -  Z-,)2 

b.  The  single  foil  interference  depends  on  the  phase  Z2) 

c.  The  multi-foil  interference  depends  on  the  phase  ( C j/  Zj  +  L,/Z,) 

The  formation  length  Zj  is  inversely  propotional  to  the  photon  energy  [Fig. 
(2.5)].  The  formation  length  Z,  of  the  foil  has  the  maximum  value  Zmax  at  E  jti  , 
(critical  energy  of  transition  radiation  above  which  the  intensity  becomes  negligible). 
E critical  can  ^*e  ^oun(^  by  taking  the  first  derivative  of  formation  length  Z„  [Eq.  (2.2)] 


F  _ - _ )[/lE 

Critical  '  j  +  y-02  >  c2  ' 


(2-17) 


where  E2  is  the  energy  corresponding  to  the  plasma  frequency  of  dielectric  foil.  For  0 
=  1/y  at  0opt  Eq.  (2.17)  becomes  (y/v^lE Zmax  can  be  found  by  substituting  Eq. 
(2.17)  into  Eq.  (2.2) 
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Figure  2.5  Formation  Length  of  Vacuum  (Ej  =  0,  y  =  1 27.  9  ^  1  y 
2hc 

Zmax  ”  E,vTr^2" 

for  0  =  l,y  at  0op[  then  Zmax  =  (v;-  E,)hcy.  For  photon  energies  E<  <  E 
formation  length  of  the  dielectric  becomes  independent  of  y  and  is 
propotional  to  the  square  of  the  medium  plasma  energy 
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Z2  =  4hc(E,  E,2)  . 
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Figure  2.6  Formation  Length  of  Al  (E,  ”  32.8  eV,  y  *  127.0—  1  y). 

At  the  critical  energy  the  formation  lengths  ,  Z,  are  approxmately  equal, 
and  the  single  interface  intensity  becomes  small.  Therefore,  above  the  critical  energy 
the  transition  radiation  intensity  will  be  negligible.  Figs.  (2.6)  and  (2.7)  show  the 
dependence  of  the  formation  length  of  Al  and  Mylar  on  photon  energy  for  0opt  —  1  y 
and  y  ■  127.  For  photon  energies  below  the  critical  energy  (2.95  keV.  1.79  keV  for  Al 
and  Mylar)  the  formation  length  is  approximately  linear  with  energy  and  independent 
of  y.  The  formation  length  Z,  of  Mylar  is  longer  than  Al  and  Mylar  has  lower  critical 
energy  than  Al  because  of  the  smaller  plasma  frequency. 

2.  Single  Interface  Intensity 

The  single  interface  intensity  d”N  dfito  [Eq.  (2.3))  limits  the  spectrum  to  the 
maximum  photon  energy  approximately  equal  to  yE,  v  2  for  a  vacuum  dielectric 
stack.  Above  this  energy,  the  spectrum  is  propotional  to  I  w4  .  Figs.  (2.S)  and  (2.9) 
show  that  d"N  dfido)  falls  off  more  quickly  for  a  Mylar  vacuum  interface  than  for  a 
Al  vacuum  interface.  The  single  interface  intensity  is  also  the  factor  which  limits 
emission  of  transition  radiation  to  a  small  forward  angle  approximately  equal  to  1  y. 
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Figure  2.7  Formation  Length  of  Mylar  (E,  *  19.9  eV,  y  =  127.  0~  1  y). 

3.  Single  Foil  Interference  Factor 

The  single  foil  interference  term  sin2(E;,Z,)  is  a  important  factor  in  the 
resonance  condition.  Single  foil  resonance  condition  occurs  at  the  photon  energy  Em 
satisfying  the  condition  Eq.  (2.10)  f ,  =  Z,  (2m- 1)71.  2  with  m  an  integer.  The  phase 
slippage  of  the  particle  during  transit  through  the  foil  is  (2m-l)7t  2.  and  the  radiation  of 
opposite  phase  produced  at  each  interface  as  the  particle  enters  and  exits  the  foil 
interferes  constructively.  The  phase  slippage  accounts  for  coherent  addition  of 
amplitudes  from  the  two  interface  of  a  single  foil  and  gives  a  peak  value  twice  as  large 
as  from  two  interface  when  the  emission  is  completly  random.  Figs.  (2.10)  and  (2.11) 
show  sinale  foil  interference  factors  of  A1  and  Mylar  with  respect  to  photon  energy . 
Mylar  has  a  longer  formation  length  Z-,  than  A1  [Figs.  (2.6)  and  (2.-)]  and  2.3  times 
bisger  foil  thickness  than  Al.  The  Mylar  has  a  second  peak  at  a  smaller  photon  energv 
than  Al.  Fig.  (2.12)  shows  the  single  foil  interference  factor  as  a  function  of  photon 
energy  for  difTerent  foil  thickness.  The  photon  energies  corresponding  to  the  single  foil 
resonance  can  be  found  from  the  first  peak  curves. 
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Figure  2.8  Single  Interface  Factor  of  A1  (y  =  127). 

4.  Muti-foil  Interference  Factor 

The  multi-foil  interference  factor  F(\1.Z)  [Eq.  (2.5)]  results  from  the  same 
interference  phenomena  as  for  the  single  foil  interference  factor.  It  includes  the  effect  of 
absorption  of  photons  within  the  foil  stack. 

As  the  number  of  foils  in  a  stack  is  increased,  intensity  at  the  fixed  photon 
energy  is  increased  up  to  the  saturated  value.  Fig.  (2.13)  shows  the  multi-foil 
dependence  according  to  the  number  of  foils.  The  absorption  factor  approximately 
increases  linearly  with  number  of  foils.  When  \I<T>  >  l  the  asymptotic  value  for 
F(M.X)  is  1  <r(<o).  After  this  value  is  attained,  the  radiation  intensity  can  not  be 
increased  by  increasing  the  number  of  foils.  A  approximate  condition  for  foil 
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Figure  2.9  Single  Interface  Factor  of  Mylar  (y  =  127). 

saturation  of  transition  radiation  spectrum  is  M<X  ^  1  .  If  the  resonance  condition  is 
satisfied,  then  the  radiation  from  successive  foils  interferes  constructively,  and  F(M.X) 
~~  M"  for  negligible  absorption  (  with  M  the  number  of  foils  ).  For  M>  >  1,  F(M.X) 
can  be  approximated  as  a  series  of  delta  functions  located  at  the  resonance  angles 
given  by  Eq.  (2.13).  Fig.  (2.14)  shows  the  photon  energy  dependence  of  F(M.X)  for 
Mylar  vacuun  with  the  number  of  foil  S  and  foil  thickness  2.5j»m  .  In  addition  to 
absorption,  particle  scattering  and  random  variation  in  foil  thickness  and  spacing 
reduce  the  resonance  effect.  (Ref.  7] 
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Figure  2.10  Single  Foil  Interference  Factor  of  A1  (6,  =  lpm). 

5.  Shape  of  the  Transition  Radiation 

The  shape  of  spectrum  [Eq.  (2.1 6>1  is  determined  by  the  three  main 
characteristics: 

a  The  absorption  of  the  radiation  in  the  foils 

b  The  critical  photon  energy  tiV)  =»  yfaj) , 

c  The  constructive  interference  between  interfaces  of  single  foil. 

The  effect  of  these  factors  is  shown  for  three  cases.  Where  the  transition 
radiation  intensity  is  plotted  as  a  function  of  energy  for  three  cases  ;  (1)  for  no 
absorption  --  single  interface  factor  [  Figs.  (2.S)  and  (2.9)  ],  (2)  for  absorption  but  no 
single  foil  coherence  -  single  interface  x  multi-foil  absorption  factor  [  Fig.  (2.15)],  (?) 


Figure  2.11  Single  Foil  Interference  Factor  of  Mylar  (£,  =  2.5fim). 

for  both  absorption  and  single  foil  coherence  --  differential  production  efficiency  [  Figs. 
(2.16).  (2.17),  (2.18).  and  (2.19)]. 


When  the  losses  due  to  the  absorption  in  the  foils  are  included,  the  reduction 
of  the  low-energy  photons  occurs  and  the  drop-off  at  the  high  energy  end  of  the 
spectrum  is  determined  by  the  critical  energy  1icoc  .  Above  this  energy  the  spectrum 
drops  off  as  (to,,  to)4  . 


The  second  term  in  Eq.  (2.16)  sin*(C-,  Z,)  accounts  for  coherent  addition  of 
amplitudes  from  the  two  interfaces  of  the  single  foil  and  gives  a  peak  value  twice  as 
large  as  from  two  interfaces  when  the  emission  is  completly  random.  This  occurs  when 
there  is  construtive  interference  between  the  waves  generated  at  the  front  and  back 
interfaces  [Ref.  8:  pp. 336-339].  The  radiation  intensity  is  maximized  when  the  thickness 
of  the  foil  is  such  that  both  the  electron  and  the  photon  travel  an  integral  number  of 
wavelengths  in  the  field  generated  at  the  first  interface. 
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Figure  2.12  Single  Foil  Interference  Dependence 
“on  Foil  Thickness  in  Mylar  (E,  =  I9.y  eV). 

The  spectral  shape  from  a  transition  radiator  is  closely  related  to  the  thickness 
of  the  individual  foil  [Eqs.  (2.4)  and  (2.5)]  not  only  because  the  formation  length  [Eq. 
(2.1)]  is  needed  for  photon  production  but  also  because  of  absorption  of  the  emitted 
radiation  in  the  foils  (Eq.  (2.6)].  For  low  energy  photons,  absorption  is  of  concern  and 
can  be  minimized  by  making  the  foils  as  thin  as  possible;  however  if  the  thicknesses  are 
made  thinner  than  the  formation  length,  photon  production  will  also  stop.  Thus,  the 
peak  production  intensity  occurs  at  the  photon  energy  corresponding  to  the  formation 
length  Z,  equivalent  to  foil  thickness  (Ref.  9:  p.3604],  There  is  an  optimum  foil 
thickness  that  balances  production  with  re-absorption  to  give  the  maximum  photon 
yield.  This  is  shown  in  Fig.  (2.20)  for  the  case  of  Mylar  3  foils.  The  number  of  foils  M 
is  kept  constant  and  the  thickness  is  varied.  The  magnitude  and  peak  emission  can  be 
predicted. 

The  spectral  shape  is  influenced  by  not  only  an  optimization  procedure  of 
minimizing  absorption  with  thin  foils  but  also  by  using  as  large  a  value  of  M  as 
feasible.  As  the  number  of  foils  increases  both  scattering  and  photon  absorption 
increase,  thereby  destroying  coherence  and  ultimately  preventing  emission  from 
photons.  Neglecting  multi-foil  scattering,  as  the  number  of  foils  increases  the  peak 
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Figure  2.13  Effect  of  the  Number  of  Foil  on  the  Multi-foil  Factor  (<y  =  0.001). 

sharpens  and  increases  in  intensity  as  M  2  [Eq.  (2.11)],  however,  since  absorption  also 
increases  with  VI  [Fig.  (2.14)]  ,  a  final  value  for  photon  production  is  reached  where 
increasing  M  does  not  result  in  increased  production.  Fig.  «'  2 . 2 1 )  shows  a  limiting 
effectiveness  of  increasing  M.  there  is  small  gain  in  intensity  after  VI  ~  20. 

The  increase  in  absorption  above  the  K-edge  results  in  a  narrower  energy 
spectrum  than  would  otherwise  be  measured.  Fig.  (2.22)  shows  the  calculated  effect  of 
K -shell  absorption  for  A1  on  the  radiation  spectrum  for  63  MeV  electron  beam  energy 
with  K-edge  1560  eV.  As  the  absorption  coefficient  increases,  the  intensity  decrease 
sharply.  The  narrow  spectrum  is  due  io  the  sudden  change  in  x-ray  absorption  at  the 
K  photoabsorption  edge  in  the  material.  The  curve  does  not  include  the  effect  of  the 
detector  resolution. 

The  purpose  of  this  experiment  has  been  divided  into  three  areas; 

(a)  The  x-ray  radiation  is  emitted  m  a  forward  cone  with  the  apex  angle  of  the 


Figure  2.14  Multi-foil  Absorption  Factor  of  Mylar  (M  **  8,  £,  =  2.5fim). 

cone  at  approximately  I  y.  The  meauremenl  of  the  apex  angle  was 
done  by  improvements  in  collecting  data  using  a  motor-driven  detector. 

The  absolute  differential  production  efficiencies  for  soft  x~rays  emitted  from 
radiator  were  compared  with  the  theoretical  pedictions  and  interpreted  in 
terms  of  absorption  and  coherence  by  measuring  the  peak  and  FIVHM 
bandwidth  energy  (full  width  of  the  spectrum  between  points  having  half  the 
maximum  production  efficiency). 

The  narrower  frequnency  spectrum  was  predicted  by  the  increase  in  absorption 
above  the  K  edge ,  and  the  effect  of  the  K  edge  was  compared  with  no  K-  shell 
absorption  edge. 
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Figure  2.15  Absorption  Factor  of  Mylar  (Single  Interface  x  Multi-foil  Absorption). 
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Figure  2.16  Calculated  Differential  Production  EfTiciencv 
from  30  Foils  Ifim  Al. 
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Figure  2  20  Calculated  Spectral  Distribution  from  8  Foils 
ot  Various  Thickness  of  Mylar  for  65  MeV. 


III.  EXPERIMENTAL  APPARATUS 
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Figure  3.1  Schematic  Diagram  of  a  Experimental  Apparatus. 

A.  EXPERIMENTAL  CONCEPT 

The  experiment  was  done  at  the  Naval  Postgraduate  School  at  Monterey 
California  using  a  medium  energy  linac.  The  experimental  apparatus  is  shown  in  Fig. 
(3.1).  Electrons  enter  from  the  left-hand  side  into  a  vacuum  chamber  where  they  pass 
through  the  foil  stacks,  and  then  through  the  dump  magnet  to  be  deflected  out  ol  the 
path  of  photon  detector.  During  their  passage,  the  x-rays  travel  entirely  in  a  10'6  Torr 
vacuum  penetrating  the  1000  A  window  of  the  proportional  counter.  Throughout  the 
run.  the  transition  radiation  photons  were  detected  with  a  gas-flow  x-ray  propotional 
counter.  Five  foil  stacks  were  mounted  on  a  movable  platform  holder  so  that  five 


different  targets  can  be  used  without  changing  the  vacuum.  A  phosphor  (ZnS)  target 
was  used  for  the  alignment  of  the  electron  beam  relative  to  the  target  chamber.  The 
detector  was  periodically  calibrated  using  a  Fe55  source.  The  current  pulse  height  was 
monitored  using  a  scintillator  just  before  the  beam  dump.  In  this  manner  the  total 
charge  was  kept  constant  for  each  measurement. 

B.  EQUIPMENT  PERFORMANCE 


Figure  3.2  Schematic  Diagram  of  a  Linear  Accelerator. 

1.  Linac 

The  linear  accelerator  (Linac)  is  used  to  produce  high-energy  electron  bunches 
[Fig.  (3.2)].  As  for  general  characteristics,  the  linac  is  similar  to  the  Stanford  linear 
accelerator  Mark  3,  which  has  1000  MeV  kinetic  electron  energy,  while  the  NPS  linac 
is  9.14  m  long  and  100  MeV.  An  electron  gun  at  the  beginning  of  the  accelerator 
injects  electrons  at  roughly  p  *  0.5.  The  design  of  the  accelerator  sections  causes  the 
energy  to  propagate  in  the  TM  mode.  The  injected  electrons  are  accelerated  to  near 
the  velocity  of  light  in  their  first  few  centimeters  of  travel.  Coupled  klystrons  produce 


a  travelling  TM  mode  wave  along  the  wave  guide.  The  electron  speed  and  the  wave 
speed  are  about  the  same,  so  that  the  electric  field  of  the  wave  will  accelerate  the 
electrons.  Thus,  the  electrons  gain  energy  at  the  expense  of  the  wave. 


Figure  3.3  Target  and  Detector  Chamber. 

2.  Target  and  Detector  Chamber 

The  target  and  detector  chamber  are  made  of  aluminum  and  the  picture  is 
shown  tn  Fig.  (3.3).  During  the  experiment  the  whole  Linac  and  a  target  chamber  are 
kept  under  vacuum  by  using  a  Turbo  Molecular  Pump  and  Mechanical  Fore  Pump. 
The  pressure  is  about  10'2 3 * * 6  Torr.  In  the  target  chamber,  there  is  a  vertically  movable 
target  ladder  on  which  the  foil  stacks  are  mounted.  The  target  ladder  position  is 
monitored  and  controlled  with  circuits  connecting  to  the  Control  Room.  The  detector 
system  will  be  discussed  in  the  next  section. 

3.  Gas  Flow  System 

The  schematic  diagram  of  the  gas  flow  system  for  the  proportional  counter  is 

shown  in  Fig.  (3.4).  Vacuum  pump  1  and  2  operate  throughout  the  experiment.  If 

valve  A  is  open  with  valves  B  and  C  closed,  the  pressure  becomes  the  same  at  valve  B. 
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Figure  3.4  Schematic  Diagram  of  a  Gas  Flow  System. 

C.  and  the  proportional  counter.  Then,  the  target  and  the  detector  chamber  are 
evacuated  simultaneously  through  the  valve  A.  Now  close  valve  A  and  open  valve  C 
and  D  with  the  vacuum  pump  1.  2  operating,  and  introduce  argon  and  propane  gas 
mixture  to  the  propotional  counter  by  opening  the  valve  B.  Set  the  valve  C  and  D  and 
valve  B  for  stable  gas  pressure  as  required.  The  pressure  meter  B  indicates  the  positive 
pressure  with  respect  to  the  target  chamber.  Next  supply  a  high  voltage  to  the 
proportional  counter  and  find  the  optimum  pressure  and  voltage  of  the  proportional 
counter. 

4.  Detection  Apparatus 

The  reason  for  the  use  of  the  proportional  counter  [Fig.  (3.5)]  as  a  detector  is 
that  the  low  photon  energy  generates  a  few  ion-electron  pairs  within  this  gas.  and  the 
generated  pulses  have  a  broad  distribution  reaching  up  to  the  amplifier  noise  level.  For 
use  in  the  soft  x-ray  region,  between  100  eV  and  3000  eV,  the  window  of  the 


Figure  3.5  Detection  Apparatus. 


proportional  counter  must  be  quite  thin  and  fragile.  The  mechanism  of  the 
proportional  counter  requires  that  the  pressure  of  the  gas  be  at  least  40  Torr.  The 
window  forms  a  built-in  filter  which  excludes  stray  electrons  and  ions  from  the  detector 
and  responses  to  the  counter  system  [Ref.  10:  p.l]. 

The  window  is  made  of  VYNS  (a  copoiymer  of  vinyl  chloride  and  vinyl 
acetate).  The  composition  is  5%  H.  41%  C,  5%  0,  49%  Cl  by  weight  with  a  5S% 
transmission  and  400  meshscreens  having  a  3x10  nun  slot  for  operating  maximum 
differential  pressure  of  S00  Torr.  For  the  propotional  counter,  the  gas  nuxture  of  90% 
argon  and  10%  propane  (P-10  gas)  is  provided  by  the  gas  flow  system  and  high  voltage 
power  (1200  V)  is  supplied.  Use  of  the  gas  molecules  as  detection  assures  l'.'O "  > 
photoionization  yield  above  the  ionization  energy,  and  with  zero  yield  below  the 
ionization  energy.  The  proportional  counter  gas  is  delivered  from  the  tank  which  will 
provides  a  choice  of  several  gas  mixtures,  allowing  a  change  in  primary  detector 
medium  by  valving  in  the  new  gas.  The  performance  depends  strongly  upon  the  gas 
pressure  which  can  be  controlled  all  the  time  during  the  experiment. 


Figure  3.6  Manson  Mode!  04  Gas  Flow  Proportional  Counter. 

The  propouonal  counter  provides  information  on  the  emitted  photon  energy 
of  radition  being  detected.  The  electronic  analysis  of  the  amplitude  of  the  output  pulse 
provides  a  technique  for  determining  the  photon  intensity.  Fig.  (3.6)  illustrates  the 
basic  design  tor  the  proportional  counter.  Since  it  is  much  easier  to  mount  a  Hat-sided 
detector  than  a  cylindrical  detector,  the  proportional  counter  is  designed  as  shown  by 
by  boring  the  cylindrical  hole  in  a  rectangular  cross  section  bar.  The  window  can  be 
mounted  by  clamping  against  the  outside  of  the  counter  body. 

Fig.  (3.6)  emphasizes  the  problem  introduced  by  this  mounting,  that  of  the 
dead  volume  inside  the  detector.  Ederer  and  Tomboulian  introduced  the  design  shown 
at  the  bottom  in  Fig.  (3.6)  for  the  use  in  the  soft  x-ray  region.  The  reason  for  the 
concern  with  the  dead  space  just  behind  the  window  is  that  the  few  electrons  released 
in  the  initial  absorption  may  be  lost  by  diffusion  back  to  the  detector  window.  The 
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Figure  3.7  Electrical  Field  Strength  and  Potential  vs.  Radius. 

counter  used  for  the  experiment  has  a  stainless  anode  50  urn  in  diameter,  maintained 
at  1000  to  1500  volts  positive  potential.  Fig.  (3.7)  shows  the  electric  field  and  voltage 
distribution  with  respect  to  the  location  from  the  anode  wire  axis.  The  anode  wire 
position  is  located  at  the  center  of  the  horizontal  scale,  which  is  a  linear  scale  of  radius 
measured  from  the  axis  of  the  counter  running  from  the  window  at  r  =  10  mm  at  the 
left  and  the  right.  The  field  increases  inversely  proportional  to  the  radios  to  very  high 
values  (94.()00  V  cm)  quite  close  to  the  anode  wire  [Ref.  10:  p. 3 ]. 

The  initial  event  ion-electron  pairs  are  separated  rapidly  by  the  field  outside  of 
the  central  region.  When  the  electrons  approach  to  within  a  few  mean  free  paths  of 
the  anode  wire,  much  energy  is  produced.  At  some  radius,  the  energy  gained  by  each 
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Figure  3.8  Proportional  Counter  Efficiency. 

event  electron  between  collisions  is  sufficient  to  cause  secondary  ionization.  Between 
the  collisions  the  energy  increases  as  the  radius  decreases  causing  most  of  the 
ionization  zone  to  be  localized  within  50  urn  of  the  anode  wire.  The  pulse  induced  on 
the  central  anode  wire  will  be  proportional  to  the  potential  difference  through  which 
the  electrons  and  ions  fall  in  their  motions  towards  the  wire  and  towards  the  grounded 
cathode  inner  surface.  The  positive  ions  move  through  most  of  the  applied  potential  in 
their  trip  to  the  grounded  cathode  and  this  positive  ion  motion  is  the  origin  of  the 
observed  output  pulse.  Since  part  of  the  trip  takes  place  in  the  very  high  field  region 
near  to  the  anode  wire  axis,  the  ion  has  a  very  high  velocity  through  the  steep  potential 
gradient  so  the  pulse  starts  off  with  a  rapid  initial  rise,  reaching  about  50°  o  of  it 
maximum  value  within  0.5fis. 

Fig.  ( 3 . S i  shows  the  counter  efficiency  with  respect  to  the  photon  energy. 
About  3.2  keV  there  is  sudden  change  in  efficiency,  but  the  energy  region  of  interest  is 
below  3.2  keV.  The  proportional  counter  can  be  moved  up  and  down  by  the  driving 
motor,  and  the  counter  always  points  at  the  virtual  center  of  the  emission  region. 

5.  Observation  Station 

The  analyzing  equipment  consists  of  an  ORTEC  142  PC  preamplifier,  an 
ORTEC  450  Research  Amplifier,  and  a  TRACER  NORTHERN  TN-7200  Pulse  Height 
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Figure  3.9  Block  Diagram  of  the  Data-coilection  Electronics. 

Analyzer  (PHA).  The  electron  beam  was  deflected  by  the  sweeping  magnet  through 
the  large  plastic  scintillator  which  served  as  the  beam  monitor.  The  PHA  received 
signal  pulses  from  a  charge-sensitive  amplifier  connected  to  the  x-ray  detector. 
Photons  were  counted  until  a  fixed  amount  of  beam  charge  had  passed  throgh  the  foil 
stack  (Fig.  (  3.91). 

The  PHA  divides  the  detection  range  into  a  predetermined  number  of 
channels,  and  then  records  graphically  the  number  of  counts  with  which  each  signal 
channel  is  detected.  The  data  is  presented  in  the  form  of  pulse  height  distribution  with 
the  number  of  counts  per  channel  displayed  on  a  linear  vertical  scale  and  the  channel 
number  displayed  on  a  linear  horizontal  scale.  The  pulse-height  is  proportional  to  the 
incident  photon  energy.  Thus  the  displayed  results  form  an  energy  dispersed  spectrum 
of  the  incident  photon  beam. 

6.  Stack  Foils  Used  during  Experiment 

Two  A1  stack-foils  and  two  Mylar  stack-foils  were  used  during  the  experiment. 
One  aluminum  stack  has  30  foils  with  I  Jim  foil  thickness,  the  other  has  1  foil  of  2f>  Jim 


Figure  3.10  Stack  Foils  Used  during  Experiment. 

toil  thickness,  one  Mylar  stack  has  S  foils  with  2.5  Jim  foil  thickness,  the  other  has  1 
foil  with  19  |im  foil  thickness.  The  Mylar  foils  were  mounted  on  ladder  2.1-mm-stick 
ladder,  while  the  A1  foils  were  mounted  on  1.5-mm-thick  ladder.  These  dimensions  were 
chosen  to  maintain  adequate  support  of  the  delicate  foils  and  keep  the  constant 
dimension  (Fig.  i  3. 10)]. 

A  stack  of  S  Mylar  foils  with  each  foil  thickness  2.5  Jim  was  constructed  [Figs. 
•  3.1  h  and  <  3.12)]  in  order  to  attain  a  coherent  spectrum.  This  stack  maintains  a  flat 
surface  with  a  spacing  1.6  mm.  Mylar  was  used  because  it  is  durable  and  does  not  tear 
easily  when  stretched. 

The  foil  stacks  produce  both  transition  and  bremsstrahlung  radiation.  To 
account  for  the  bremsstrahlung  background  for  each  radiator,  a  single  foil  with 


Figure  3.1 1  Design  of  Mylar  Coherent  Stack- 


thickness  equivalent  to  stack  foils  was  constructed.  The  measured  bremsstrahlung 
generated  from  the  single  foil  using  equal  beam  current  was  then  subtracted  from  that 
produced  by  the  corresponding  foil  stacks. 


IV.  RESULTS  AND  CONCLUSIONS 


A.  RESULTS 

Measurements  of  the  total  number  of  counts  were  carried  out  for  several  different 
positions  of  the  detector.  All  measurements  were  for  the  same  counting  interval.  The 
detector  plane  was  49"  from  the  stack  target.  Motion  of  the  detector  was  monitored 
by  a  counter  with  a  scale  of  0.0124"  digit.  The  detector  was  on  the  beam  axis  at  a  scale 
reading  of  335.  Table  1  shows  the  total  number  of  counts  in  region  of  interest  for 
several  different  detector  positions  for  65  Mev  electrons  incident  on  30  A1  stack  foils. 

TABLE  1 

MEASUREMENTS  OF  TOTAL  COUNTS 


Dial  No. 

Distance 

Cone  Angle 

1 

Counts  No. 

285 

0.620" 

12.65mr 

10102 

! 

290 

0.558" 

1 1.39mr 

10112 

295 

0.496" 

10.12mr 

10123 

300 

0.434" 

8.86mr 

10134 

i 

305 

0.372" 

7.59mr 

i 

10144 

310 

0.310" 

6.33mr 

i 

10185 

315 

1 

0.24S"  j  5.<>6mr 

10165 

| 

45 


TABLE  2 
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Figure  4.6  Measured  Pulse  Height  Counts 
from  1  Foil  I9jim  Mvlar. 


the  measured  foil  stack  spectrum.  The  theoretical  plots  were  calculated  by  using  Eq. 
(2.16). 

The  ratio  of  the  transition  radiation  to  that  of  bremsstrahlung  was  2.3  for  A1  and 
2.7  for  the  Mylar  foil.  These  ratios  are  lower  than  expected  from  the  theory  and  the 
reason  may  be  because  there  are  other  sources  of  background  which  were  originally 
attributed  to  bremsstrahlung.  A  sudden  change  in  intensity  was  observed  at  energy  3.1 
keV  for  Mylar  and  it  is  due  to  the  counter  efficiency  change  from  0.2  to  0.5  [Fig.  (3.8)]. 

Table  3  is  a  listing  of  four  foil  stacks,  their  respective  number  of  foils,  calculated 
peak  frequencies,  measured  frequencies,  calculated  FVVHM,  and  measured  FVVHM. 
For  three  of  the  foil  stacks,  the  measured  photon  peak  positions  are  found  to  agree 
with  the  calculated  values.  The  single  foil  Mylar  is  noticeably  different  because  of  big 
fluctuations  in  small  photon  energy  and  of  the  modelling  of  unknown  background. 
Comparison  between  the  theoretical  curves  and  their  respective  experimental 
measurements  are  limited  to  the  position  of  the  peak  and  the  bandwidth.  The 
measured  FWHM  energy  for  A1  stack  is  reduced  by  150  eV,  and  it  is  due  to  bandwidth 
narrowing  associated  with  K-edge  absorption  (E=  1.54  keV). 

B.  CONCLUSIONS 

The  absolute  differential  production  efficienies  (photons/eV  sr  electron  )  for  soft 
x-rays  emitted  from  each  of  four  transition  radiator  were  measured  for  incident 
electron-energy  of  65  MeV.  The  radiators  were  made  of  stacks  of  1.0-pm-thick  30  foils 
Al.  25.4-ftm-thick  single  foil  Al,  2.5-jim-thick  8  foils  Mylar,  and  19-(im-thick  single  foil 
Mylar.  The  radiation  spetra  were  most  intense  between  0.6  and  1.8  keV  for  Als.  0.9  and 
3.1  keV  for  Mylar,  peaking  at  1.10,  1.00,  1.44,  and  1.58  keV  respectively. 

The  transition  radiation  cone  was  measured  in  different  ways,  and  the 
measurement  is  in  agreement  with  optimum  apex  angle.  The  radiation  peak  can  be 
selected  by  proper  selection  of  foil  number,  and  thickness.  For  the  case  of  foils 
separated  by  a  vacuum  ,  the  peak  frequency  and  bandwidth  are  dependent  upon  the 
foil  thickness  and  plasma  frequency.  Measured  peak  energies  correspond  to  the  theory. 
Aluminum  bandwith  narrowing  associated  with  K-edge  absorption  was  observed. 

The  highly  directional  beam,  high  photon  production  make  an  elecrtron  beam 
competitive  as  an  x-ray  source,  especially  for  applications  using  a  low  energy,  low 
current  electron  beam  and  requiring  a  small  x-ray  beam  size.  X-ray  lithography,  x-ray 
microscopy,  and  other  medical  application  are  possible. 


50 


However  refinements  remain  to  be  done.  Some  specific  suggestions  are  given 
below. 

1.  The  magnet  was  used  to  deflect  electrons  from  the  photon  path,  but  a  part  of 
them  entered  the  photon  path.  They  made  noise.  More  powerful  magnets  are 
needed  to  eliminate  electrons. 

2.  The  experimental  area  was  very  noisy.  It  affected  the  observed  spectrum.  Noise 
could  be  decreased  by  electrically  shielding  the  detector,  which  would  decrease 
the  electric  field  noise.  Magnetic  field  noise  associated  with  the  LINAC 
klystrons  could  be  reduced  by  thick  concrete  shielding. 

3.  The  pipe  through  which  photons  pass  is  about  1.25",  and  the  detector  has  the 
dimension  l/8"x3/8"  of  slotted  window.  The  resolution  smearing  occured  by  the 
shadow  effect  of  photons.  An  x-ray  detector  with  very  fine  resolution  would  be 
required. 

Continued  research  of  transition  radiation  include  the  following  for  future  work: 

1.  Investigation  of  the  effect  of  multi-foil  scattering  in  transition  radiation 

2.  Determination  of  the  spectrum  distribution  of  a  radiator  having  the  several 
different  periodic  media  (such  as  AI,  air.  and  Mylar) 

3.  Determine  the  dependency  of  the  intensity  and  frequency  spectrum  on  the 
resonance  peak 


APPENDIX  A 

ELECTRON  LINEAR  ACCELERATOR  PARAMETERS 


TABLE  4 

LINAC  PARAMETER 


Max.  Enerav 

120  MeV 

Max.  Average  Current 

20  jiamps 

Length 

9.14  m 

Pulse  Repetition  Frequency 

60  Hz 

Pulse  Width 

1  ns 

Operation  Frequency 

2.356  GHz 

Pulse  Repetition  Rate 

120  PPS 

Pulse  Duration 

2.5  fis 

j 

Nr.  of  Klystron  j  3 

1 

Peak  Power  per  Klystron 


21  MW 


APPENDIX  B 

CALCULATION  OF  PLASMA  FREQUENCY 


TABLE  5 

PLASMA  FREQUENCY 


material 

CO 

p 

A1 

32.8  eV 

Be 

•  24.5  eV 

Cu 

54.6  eV 

Mylar 

19.9  eV 

_ 1 

In  the  absence  of  a  magnetic  field,  a  plasma  rasonates  electrostatically  with  the 
frequency.  Ignoring  the  ion  motion,  the  plasma  frequency  is  given  by 


to  =  Vne2  cr 


m 


(2.1) 


where  e  is  1.6xl0'!9  (coul),  h  is  1.055xl0’34  (J.sec),  m  is  9. 1  x  10*3 1  (kg),  and  cQ  « 
8.85xl0‘12  (F/m)  .  Substituting  these  values  into  Eq.  (2.1),  the  plasma  frequency 
becomes 


hcop  =  ^/n"  3.7I8xIO_J  1  (eV) 


(2.2) 


The  electron  concentration  can  be  obtained  bv 


where  n  is  the  number  of  electrons  per  cm3,  Z  is  atomic  number,  N  is  avogadro 
number,  p  is  density  (gm/cm3),  and  A  is  atomic  weight.  For  A1  n  =  7.8xl023 
(electrons; cm3),  the  plasma  frequency  has  the  value  tko  =  32.84  (eV). 
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APPENDIX  C 


MASS  ATTENUATION  COEFFICIENTS  FOR  PHOTONS  IN  AL  AND 

MYLAR 


Figure  C.l  Mass  Attenuation  Coefficient  for  Photons  in  Al. 

Corresponding  linear  coefficients  for  Al  can  be  obtained  by  multiplying  all  curves 
by  p  -  2.70  (g  em3)  for  Al  [Fig.  (C.l)|  and  p  *  1.395  (g  cm3)  for  Mylar  [Fig.  <C.2i| 
(at  f  *  0°C,  P  ■  760  mmHg).  Mylar  (CjH402)  mass  attenuation  coefficients  are 
obtained  approximately  by  taking  p  p  -  1  3(|»  p)0  •*-  5  3( ft  p)  c. 


Figure  C.2  Mass  Attenuation  Coefficient  for  Photons  in  Mylar. 
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APPENDIX  D 

BASIC  PROGRAM  (HP  COMPUTER) 


INFORMATION  . . . . . . . 

THE  OBJECTIVE  OF  THIS  PROGRAM  IS  TO  CALCULATE  ANO  ORAM 
THE  ABSOLUTE  DIFFERENTIAL  PRODUCTION  EFFICIENCIES  FROM 
STACK  MATERIAL.  THE  ORAHING  METHOO  IS  THE  FUNCTION 
CURVE  OF  DISSPLA. 

hhhhhhhhhhhw  VARIABLE  DEFINITION  . . mhhhhhhi 

L  :  FOIL  THICKNESS  IN  Ua  • 

X  :  ENERGY  IN  KEV 

K  :  ENERGY  IN  FREQUENCY 

0  :  ABSOLUTE  COEFFICIENT 

Y  :  MULTI-FOIL  ABSORPTION  COEFFICIENT 

0  :  PLASMA  FREQUENCY  IN  KEV 

A  :  PLASMA  FREQUENCY  IN  FREQUENCY 

N,AA,B  :  SINGLE  INTERFACE  FACTOR 

U  :  FORMATION  LENGTH  Z1 

G,E  :  FORMATION  LENGTH  Z2 

F  :  SINGLE  FOIL  FACTOR 

GG>H  :  TOTAL  PRODUCTION  EFFICIENCY 

M  :  TOTAL  ABSORPTION  FACTOR 

mwi»wimin>»i»  VARIABLE  DEFINATION  ihiwhhhhhhmhmhmmmmhhwhhmhhhmhhhi 


A  PROGRAM  FOR  PLOTTING  BY  SHERPAOR  BY  TEK618 

C  GRAPH  OF  XY  DATA  GIVEN 

PROGRAM  GRFZO 
CALL  TEK6I8 

C  CALL  PRTPLT 

CALL  XYPLOT 
CALL  DONEPL 
STOP 
ENO 

C  PLOTTING  THE  GIVEN  FUNCTION 
SUBROUTINE  XYPLOT 


MAIN  PROGRAM 


DIMENSION  X(  300  )  ,Y(  300  1 ,01  300  I  ,K<  300  )  ,At  300  )  ,B(  300  )  ,C(  300  )  ,E(  300  ) 
DIMENSION  F<  300  I  ,G(  300  I  >H(  300  )  .Ml  300  )  .LI  300  )  >NI  300  )  .01  300  > .PI  300  ) 
DIMENSION  Qt 300  I ,R(  300  I .SI  300 ) ,T< 300 ) ,U( 300 ) .VI 300 ) ,H( 300  )  .AAI 300  > 
DIMENSION  BBI 300  I  .CCI 300  I  ,FF<  300  ) , EE  I  300  I  >GGI 300  )  ,YY( 300  I 
REAL  X.Y.O.K.A.B.C  .E  .G.H.M.L.Z.N.O.P.Q.R.S.T .U.V.H.AA.BB.CC.FF >GG 
C  REAL  YY 

INTEGER  I.J 
CALL  PAGE! 12. ,9.  I 
CALL  AREA20I 9. ,5.9 1 
CALL  AREA20I 5.5.4.  ) 

CALL  YAXANOI 90.  ) 

CALL  XASANGIO.  I 

CALL  MX2ALFI  ,L/CSTD,.V'  ) 

CALL  »<3ALFI 'ITALIC', 'Z'  ) 

CALL  HEIGHT! 0.2) 

CALL  XNAMEI  ‘PHOTON  ENERGY  I KEV) ',20) 

CALL  XNAMEI 'NUMBER  OF  FOIL'  ,19) 

CALL  YNAME 1 'PHOTONS  /  EV/STERAD . /ELECTRON )  10**< -9  )  • ,39  ) 

CALL  YNAME I 'SINGLE  INTERFACE  *  MULTI-FOIL  AB.  10**l -8  ) ' ,92  ) 

CALL  YNAME I 'ABSORPTION  COEFFICIENT  I  CM**2  /GM )  IN  LOG', 91) 

CALL  YNAME! 'MULTI-FOIL  INTERFERENCE  FACTOR', 30) 

CALL  COMP LX 

CALL  HEAOINI 'PRODUCTION  EFFICIENCY ' ,21 ,2 . 5,1 ) 

CALL  HEAOINI 'MULTI-FOIL  INTERFERENCE ' ,23,1.5,1 ) 

CALL  RESETI 'HEIGHT'  ) 

CALL  FRAME 
CALL  CROSS 

CALL  GRAFIO. ,100. ,900. ,0. ,100. ,500. I 

CALL  DOT 

CALL  GRID)  1,1) 

CALL  RESETI 'ALL'  ) 

CALL  POLYS 
CALL  THKCRVI .05) 

CALL  MARKER!  1) 


REPRODUCED  AT  GOVERN ME U  T  F^rru'.r 


CALL  HEIGHT!.  15) 

Z  =  0. 

00  30  J  *1,5 
Z  *  Z  ♦  1 
DO  20  I  >  1.300 
M!  I  )  *  4*Z 
L(  X 1  *  1 

XII)  *  FLOAT!  1 1/25 
Kl  I )  =  1.51fc59*XI I  ) 

Oil)*  3.91-0. 39*X(  I  ) 

01  I  )  *  5.55-XI  I  1*0.55 
01  I  )  *LOG10I  XI I  I  I 
£EI  I  >  *-2.1*LOG10IX(III*2.S 
FFII)  =  -3.0792*L0G10(XID  1*5.7024 
Dll)  *  10**(  E£l  I  I  1*0. 003*2 . 7 
EEIII  =  -2.1*LOG10IX(II  1*2.5 
Oil)  *  10**1  EEl  I )  l»l*l£-4«2 . 7 
01 1  I  *  0.001 

VI  X  I  *  1/01  I  1*1  1-EXPI -XII  1*011 ) )) 

0(1  )  *  30.9E-S 
All)  *  1 .5165°*0(  I  ) 

BID  *  1 .8282E-7*Kt  I  )**(  -5  )*A!  I  )**4 
C(I)  *  1 . 2S*E  -4*1  A 1 1  l/KI  I  )  I**  2 
Ell)*  l/KI  I  1*1 . 2E 9/C!  I  ) 

UII  )  *  6.4E-6/XI I ) 

HID  *  (U! I  l-EI  1  I*1E-18  l**2 

AAlD*  3 . 174E-26*HI  I  )*KI  I  I»1E18*1 . 5172E1S*1E4*4 
Fill  «  SINI  II I  l*lE-6/E<  I  I*1E18  1**2 
GII)  *  80«0.6/l  1 . 236E-4*I  Al  1  l/KI  I  I  l**2  ) 

Hill  S  GII  l»»2»F(I  l*B(II*¥II  I*1E  -18*1 . 523E15 
Mil)  *  YII  l»F  1 1  I 

PID  *  1.23bE-4*(1.236E-4*(OIII/X(II  1**21 

Ni  I)  *  1 ,8282£-25*0l  I  l**4»X(  I  1**1  -5  )*P(  I  l**l  -2  )*1E17 

Til)  *  S175/UI I  l*LI  I  i/EI  I  )«IE12 

0(11  *  I  Slr!l8»T(  I  I  /SINI  Tl  I  I  I  l««2 

Til)  *  Li  1  l/EI  I  I*1E12 

0(11  *  1*EXP( -8*01  II  >-2«EXPl -4*0(1  >  >*C0S(  16*TI  I )  > 
RID  *  1*EXPI -01 1)  l-2*EXPI -01 1  1/2  l*C0S(  2*T(  II  I 
SI  II  *01  I  l/RII  ) 

GGI  I  I  *AA(  I  l*Y 1 1  l*FI  I ) 

CCI  I  I*  AAI  I  l*SI  I  I 

BBI  I  I*  AAI  I  l*FI  I  l*YI  I  I 

VI II  *  NI  I  l*F(  I  l*St  I  )*1.523E15 

WRITE  16.55  IXII  I.GGII) 

55  FORMAT  I  •  * .4E10.3  ) 

REAOI  b  .55  IX.Y1.Y2 
55  FORMAT  1011. 3. 2X. 2012. 3) 

20  CONTINUE 

CALL  CURVE! X.Q, 300,0 ) 

CALL  EN0GRI0) 

30  CONTINUE 

CALL  RESET  I  'ALL*  I 
CALL  EN0PLI0I 
CALL  ENDGRIO) 

RETURN 

ENO 


5S 


'j 
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